Search for: All records

Abstract We present nucleosynthesis and light-curve predictions for a new site of the rapid neutron capture process (r-process) from magnetar giant flares (GFs). Motivated by observations indicating baryon ejecta from GFs, J. Cehula et al. proposed that mass ejection occurs after a shock is driven into the magnetar crust during the GF. We confirm using nuclear reaction network calculations that these ejecta synthesize moderate yields of third-peakr-process nuclei and more substantial yields of lighterr-nuclei, while leaving a sizable abundance of free neutrons in the outermost fastest expanding ejecta layers. The finalr-process mass fraction and distribution are sensitive to the relative efficiencies ofα-capture andn-capture freeze-outs. We use our nucleosynthesis output in a semianalytic model to predict the light curves of novae breves, the transients following GFs powered by radioactive decay. For a baryonic ejecta mass similar to that inferred of the 2004 Galactic GF from SGR 1806-20, we predict a peak UV/optical luminosity of ∼10³⁹–10⁴⁰erg s⁻¹at ∼10–15 minutes, rendering such events potentially detectable to several Mpc following a gamma-ray trigger by wide-field transient monitors such as ULTRASAT/UVEX. The peak luminosity and timescale of the transient increase with the GF strength due to the larger ejecta mass. Although GFs likely contribute 1%–10% of the total Galacticr-process budget, their short delay-times relative to star formation make them an attractive source to enrich the earliest generations of stars. 

Abstract The origin of heavy elements synthesized through the rapid neutron capture process (r-process) has been an enduring mystery for over half a century. J. Cehula et al. recently showed that magnetar giant flares, among the brightest transients ever observed, can shock heat and eject neutron star crustal material at high velocity, achieving the requisite conditions for anr-process. A. Patel et al. confirmed anr-process in these ejecta using detailed nucleosynthesis calculations. Radioactive decay of the freshly synthesized nuclei releases a forest of gamma-ray lines, Doppler broadened by the high ejecta velocitiesv ≳ 0.1cinto a quasi-continuous spectrum peaking around 1 MeV. Here, we show that the predicted emission properties (light curve, fluence, and spectrum) match a previously unexplained hard gamma-ray signal seen in the aftermath of the famous 2004 December giant flare from the magnetar SGR 1806–20. This MeV emission component, rising to peak around 10 minutes after the initial spike before decaying away over the next few hours, is direct observational evidence for the synthesis of ∼10⁻⁶M_⊙ofr-process elements. The discovery of magnetar giant flares as confirmedr-process sites, contributing at least ∼1%–10% of the total Galactic abundances, has implications for the Galactic chemical evolution, especially at the earliest epochs probed by low-metallicity stars. It also implicates magnetars as potentially dominant sources of heavy cosmic rays. Characterization of ther-process emission from giant flares by resolving decay line features offers a compelling science case for NASA’s forthcoming COSI nuclear spectrometer, as well as next-generation MeV telescope missions. 

Abstract Some hydrogen-rich core-collapse supernovae (SNeIIP) exhibit evidence of a sustained energy source powering their light curves, resulting in a brighter and/or longer-lasting hydrogen recombination plateau phase. We present a semi-analytic SNIIP light-curve model that accounts for the effects of an arbitrary internal heating source, considering as special cases⁵⁶Ni/⁵⁶Co decay, a central engine (magnetar or accreting compact object), and shock interaction with a dense circumstellar disk. While a sustained internal power source can boost the plateau luminosity commensurate with the magnitude of the power, the duration of the recombination plateau can typically be increased by at most a factor of ∼2–3 compared to the zero-heating case. For a given ejecta mass and initial kinetic energy, the longest plateau duration is achieved for a constant heating rate at the highest magnitude that does not appreciably accelerate the ejecta. This finding has implications for the minimum ejecta mass required to explain particularly long-lasting SNe, such as iPTF14hls, and for confidently identifying rare explosions of the most massive hydrogen-rich (e.g., Population III) stars. We present a number of analytic estimates that elucidate the key features of the detailed model. 

Abstract Relativistic jets from a Kerr black hole (BH) following the core collapse of a massive star (“collapsar”) is a leading model for gamma-ray bursts (GRBs). However, the two key ingredients for a Blandford–Znajek-powered jet—rapid rotation and a strong magnetic field—seem mutually exclusive. Strong fields in the progenitor star’s core transport angular momentum outward more quickly, slowing down the core before collapse. Through innovative multidisciplinary modeling, we first use MESA stellar evolution models followed to core collapse to explicitly show that the small length scale of the instabilities—likely responsible for angular momentum transport in the core (e.g., Tayler–Spruit)—results in a lownetmagnetic flux fed to the BH horizon, far too small to power GRB jets. Instead, we propose a novel scenario in which collapsar BHs acquire their magnetic “hair” from their progenitor proto–neutron star (PNS), which is likely highly magnetized from an internal dynamo. We evaluate the conditions for the BH accretion disk to pin the PNS magnetosphere to its horizon immediately after the collapse. Our results show that the PNS spin-down energy released before collapse matches the kinetic energy of Type Ic-BL supernovae, while the nascent BH’s spin and magnetic flux produce jets consistent with observed GRB characteristics. We map our MESA models to 3D general-relativistic magnetohydrodynamic simulations and confirm that accretion disks confine the strong magnetic flux initiated near a rotating BH, enabling the launch of successful GRB jets, whereas a slower-spinning BH or one without a disk fails to do so. 

Abstract UsingAthena++, we perform 3D radiation-hydrodynamic calculations of the radiative breakout of the shock wave in the outer envelope of a red supergiant (RSG) that has suffered core collapse and will become a Type IIP supernova. The intrinsically 3D structure of the fully convective RSG envelope yields key differences in the brightness and duration of the shock breakout (SBO) from that predicted in a 1D stellar model. First, the lower-density “halo” of material outside of the traditional photosphere in 3D models leads to a shock breakout at lower densities than 1D models. This would prolong the duration of the shock breakout flash at any given location on the surface to ≈1–2 hr. However, we find that the even larger impact is the intrinsically 3D effect associated with large-scale fluctuations in density that cause the shock to break out at different radii at different times. This substantially prolongs the SBO duration to ≈3–6 hr and implies a diversity of radiative temperatures, as different patches across the stellar surface are at different stages of their radiative breakout and cooling at any given time. These predicted durations are in better agreement with existing observations of SBO. The longer durations lower the predicted luminosities by a factor of 3–10 (L_bol∼ 10⁴⁴erg s⁻¹), and we derive the new scalings of brightness and duration with explosion energies and stellar properties. These intrinsically 3D properties eliminate the possibility of using observed rise times to measure the stellar radius via light-travel time effects. 

Abstract We explore the three-dimensional properties of convective, luminous (L≈ 10^4.5–10⁵L_⊙), hydrogen-rich envelopes of red supergiants (RSGs) based on radiation hydrodynamic simulations in spherical geometry usingAthena++. These computations comprise ≈30% of the stellar volume, include gas and radiation pressure, and self-consistently track the gravitational potential for the outer ≈3M_⊙of the simulatedM≈ 15M_⊙stars. This work reveals a radius,R_corr, around which the nature of the convection changes. Forr>R_corr, though still optically thick, diffusion of photons dominates the energy transport. Such a regime is well studied in less luminous stars, but in RSGs, the near- (or above-)Eddington luminosity (due to opacity enhancements at ionization transitions) leads to the unusual outcome of denser regions moving outward rather than inward. This region of the star also has a large amount of turbulent pressure, yielding a density structure much more extended than 1D stellar evolution predicts. This “halo” of material will impact predictions for both shock breakout and early lightcurves of Type IIP supernovae. Inside ofR_corr, we find a nearly flat entropy profile as expected in the efficient regime of mixing-length theory (MLT). Radiation pressure provides ≈1/3 of the support against gravity in this region. Our comparisons to MLT suggest a mixing length ofα= 3–4, consistent with the sizes of convective plumes seen in the simulations. The temporal variability of these 3D models is mostly on the timescale of the convective plume lifetimes (≈300 days), with amplitudes consistent with those observed photometrically. 

Abstract We present the discovery of the Type II supernova SN 2023ixf in M101 and follow-up photometric and spectroscopic observations, respectively, in the first month and week of its evolution. Our discovery was made within a day of estimated first light, and the following light curve is characterized by a rapid rise (≈5 days) to a luminous peak (M_V≈ − 18.2 mag) and plateau (M_V≈ − 17.6 mag) extending to 30 days with a fast decline rate of ≈0.03 mag day⁻¹. During the rising phase,U−Vcolor shows blueward evolution, followed by redward evolution in the plateau phase. Prominent flash features of hydrogen, helium, carbon, and nitrogen dominate the spectra up to ≈5 days after first light, with a transition to a higher ionization state in the first ≈2 days. Both theU−Vcolor and flash ionization states suggest a rise in the temperature, indicative of a delayed shock breakout inside dense circumstellar material (CSM). From the timescales of CSM interaction, we estimate its compact radial extent of ∼(3–7) × 10¹⁴cm. We then construct numerical light-curve models based on both continuous and eruptive mass-loss scenarios shortly before explosion. For the continuous mass-loss scenario, we infer a range of mass-loss history with 0.1–1.0M_⊙yr⁻¹in the final 2−1 yr before explosion, with a potentially decreasing mass loss of 0.01–0.1M_⊙yr⁻¹in ∼0.7–0.4 yr toward the explosion. For the eruptive mass-loss scenario, we favor eruptions releasing 0.3–1M_⊙of the envelope at about a year before explosion, which result in CSM with mass and extent similar to the continuous scenario. We discuss the implications of the available multiwavelength constraints obtained thus far on the progenitor candidate and SN 2023ixf to our variable CSM models.